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We try to apply a constituent quark model (a variety chiral constituent quark model) and the 
resonating group approach for the multi-quark problems to compute the effective potential between 
the NN in S'-wave (the quarks in the nucleons N and N, and the two nucleons relatively as well, are 
in S wave) so as to see the possibility if there may be a tight bound state of six quarks as indicated 
by a strong enhancement at threshold of pp in J/tp and B decays. The effective potential which we 
obtain in terms of the model and approach shows if the experimental enhancement is really caused 
by a tight 5*-wave bound state of six quarks, then the quantum number of the bound state is very 
likely to be I = f , J PC = CT+. 
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I. INTRODUCTION 

Recently the BES collaboration in the radiative decay 
J/tp — ► "/pp observed a sharp enhancement at threshold 
in the pp invariant mass spectrum[lj . They, further based 
on a multi-quark bound state conjecture, tried to fit the 
enhancement by means of an S-wave Breit-Wigner res- 
onance function, and obtained the resultant mass peak 
at M = 1859tf (stat)tf 5 (sys) MeV and width T less 
than 30 MeV. In the meanwhile, Belle also reported their 
observations that in the decays B + — > K + pp and 
an enhancement in the pp invariant mass 
distribution near the threshold was also seen clearly. 

Theoretical investigations of baryon-antibaryon bound 
states date back to the proposal of Fermi and Yang0| to 
make the pion with a nucleon-antinuclcon pair. The ex- 
tensive and excellent reviews are given in Ref [j| . In tra- 
ditional nuclear interaction theories within the potential 
framework which are based on the single meson exchange 
mainly, it is shown that iViV-system is more attractive 
than AW-system due to the fact that in the theories there 
is a strong w-exchange. Therefore possible bound states 
or resonances of a nucleon-antinucleon system (the so- 
called nuclear baryonia) have been proposed for many 
years. 

Recently being encouraged by the observations and 
with intuition, several interpretations on the observa- 
tions near and below the baryon-antibaryon threshold 
were suggested: (i) Based on the fact that the radiative 
J/ijj decay may be a gluon-rich process, L. RosnerQ ex- 
plained the enhancement is due to the iso-scalar state 
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with J PC = ±+ being coupled to a pair of gluons; (ii) 
B.S. Zou et al with a K-matrix approach 7] summed up 
the one-pion-exchange final state interaction and showed 
an important contribution to make the enhancement be- 
havior near the threshold; (iii) Many phenomenological 
models were used to explain this anomalous (sharp) en- 
hancement as a bound state or a resonance, e.g., A. Datta 
et al @ accounted it for a new bound state in a simple 
potential model with a A • A confining interactions; X.-G. 
He et al@ also gave their explanation by using linear a 
model. 

The traditional interaction theories to study NN sys- 
tem, such as boundary condition model, absorptive po- 
tential model, optical model and coupled-channel models 
etc, specially put forward their own method to handle 
the short-range part of NN, although in comparatively 
long range t hey are quite similar. For instance, the op- 
tical model 0, 0, [3 E3 provides a realistic picture 
of the scattering process by introducing an imaginary 
part of the potential to reproduce the effect of other 
excitation channels, and the spin-isospin dependence of 
meson-exchange and channel dependence in the annihi- 
lation potential as well are also taken into account. The 
coupled-channel model either describes annihilation in 
terms of baryon exchange with the same baryon-meson 
coupling as in the Yukawa potential |14| or uses semi- 
phenomenological potential adding partial-wave-analysis 
to study NN interaction[H[ll[l3. Although the agree- 
ment with collision experiments is obtained, from the 
QCD point of view each of them has shortcoming re- 
spectively. For example, it is hard to imagine that a 
baryon-exchange picture can be applied to such a short- 
range where quarks are 'overlapped' and the color octet 
configuration must be considered. 

In fact, the early studies within the traditional me- 
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son exchange (mesons exchange between the nucleons as 
whole) framework found that, if neglecting annihilation 
channels i.e. taking the real part of the effective po- 
tential only, many bound states might be formed, while 
annihilation effects i.e. the imaginary part of the poten- 
tial, were included, the binding force decreased and some 
bound states were washed out [l^. Moreover, in the 
earlier days experimentally the data there was no clear 
evidence to imply the existence of strongly bound states. 
Therefore, we suspect that the bound state should be 
tight i.e. the two components have great overlapping, if 
a pair of NN really form a bound state as indicated by 
BES and Belle experiments, i.e., the picture with meson- 
exchanges between the 'whole' nucleons to describe the 
interaction may not contain all of the key effects, at least, 
the behavior of TV TV interaction in short-range should be 
re-considered carefully with a novel way. 

A possible way out for the tight bound state prob- 
lem is Quantum chromodynamics (QCD), the fundamen- 
tal strong interaction physics. Namely we should start 
with the quark level, rather than that of the baryon 
level. QCD has been already proved to be the right ap- 
proach at high energy. At low energy, because of the 
non-perturbative nature of QCD, one has to rely on ef- 
fective theories underlining QCD and/or QCD-inspired 
models to get some insight into the phenomena of the 
hadronic world. The constituent quark model is one of 
them. It has achieved a lot of successes in describing 
spectrum for single baryon, the baryon-baryon interac- 
tions and the binding of two baryons such as the nucleus 
deuteron, even having found possible strong binding of 
fTfr systems [H Hfl, HI El H3| etc. Therefore, ex- 
panding the constituent quark model into NN system is 
an interesting tentative. 

In fact, in order to fitting all of the data from single 
baryon spectrum to interactions between two baryons, in 
literature the constituent quark model has variety ver- 
sions. In this paper, to study of the possibility of NN 
bound states as indicated in Jji\) — > jpp, B + — > K + pp 
and B° — ► D°pp decays within the framework of resonat- 
ing group method, we take a variety chiral quark model, 
which may fit the data of single baryon and deuteron 
etc. Besides including n, a and one-gluon exchange, in 
the concerned model, we take into account the contribu- 
tions from the annihilations of a pair quarks into a meson 
accordingly and into one gluon as well. In order to keep 
the model well-described for baryon spectrum and NN 
scattering data, all of the model parameters have been 
fixed as possible as those related to NN interaction and 
baryon spectrum etc. 

The paper is organized as follows. In Section II, we 
derived the formalism, then gave the Hamiltonian and 
wave function, after that, gave a brief summary about 
resonating group method, which is used to do dynamical 
calculation. In Section III, we show the results and gave 
discussion. 



II. HAMILTONIAN, WAVE FUNCTION AND 
RESONATING GROUP METHOD 

A. Formalism 

We think the basis of the constituent quark model for 
light flavor systems would be better to set on effective 
chiral quark model (with necessary extensions) of QCD. 
Namely the dynamics should be described by the effective 
chiral quark Lagrangian (pseud-goldstone particles and 
light quarks as active degrees of freedom in the theory) 
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Since at present we consider the NN system, so for the 
Lagrangien we restrict ourselves to consider two flavors 
only. As mentioned in Introduction, in order to fit nu- 
clear data some 'modification' on the quark chiral model 
is needed, such as to add a scalar a 'meson' [25| and 
a 'survived' effective gluon into the model etc. Thus 
with the modification the relevant quark-gluon-meson 
Lagrangian can be re-written as 

L c h = $q(i ft - m q )^ q + g c htpq{<7 + «75^ ■ ~^)ip q , (6) 

where ip q denotes quark field, and m q (q = u, d) are the 
constituent quark masses. Current masses turns to con- 
stituent masses accordingly is a consequence of the chiral 
symmetry breaking. g c h is the vertex coupling (quarks to 
mesons) constant and D M — + ig s ^-A afl (a is an in- 
dex of color space and A a is Gell-Mann matrix for SU C (3) 
color) A atJ - is gluon field. To study the baryon structure 
and baryon-baryon interaction, the survived gluon field 
jyifi j s introduced so as to take care of the necessary non- 
perturbative contributions in the model. In addition, to 
provide the non-perturbative QCD effects at long dis- 
tance, an effective confinement in potential is needed. 
We will discuss them later. 

Here when computing the effective potential between 
N and N, we only take into account the tree Feynman di- 
agrams: the relevant exchange and annihilation diagrams 
(see Fig. 1-3). 
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FIG. 1: The exchange diagrams between two quarks. 
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FIG. 2: The exchange diagrams between quark and anti- 
quark. 



The effective potential derived from the 'exchange' be- 
tween quarks has been considered by many authors such 

as those of Refs. is us in mm, now we should con- 
sider those derived from the exchange and annihilation 
between quark and antiquark. 

(i). The effective potential derived from 'exchange': 
According Feynman rule, the T-matrix of cr-exchange 
between quark and anti-quark can be written as 



u(p'i,s[)u(pi,si)- 



9lhv{P2, s 2 )v{p' 2 , s' 2 ) . (7) 



Here, we have inserted a form factor F(l? ) to cut off the 
contributions from the short-distance region where the 
chiral symmetry has not broken. Generally, it takes the 



form F(lf ) = ( 



A 2 +T 



a/2 



Pi, Si,p[, s[ are the four- 



vector momenta and spin z-projections of initial quark 
and final quark, respectively. £>2, S2,P2' s 2 are the four- 
vector momenta and spin z-projections of initial anti- 
quark and final antiquark, respectively, u and v are as- 
sumed to be free Dirac spinors for quarks and antiquarks, 
respectively, but soon in the paper we will replace them 
by two component Pauli spinors for convenience because 
of the non-relativistic nature for the problem. 

When we take the non-relativistic limit to turn the 
Dirac spinors to Pauli ones, there is an additional mi- 
nus sign from the interchange of the effective interaction 
uuuu to the uuvv so it cancels the extra fermion minus 
sign. Therefore the cr-exchange potential is universally 
attractive, no matter between a pair of fermions or be- 
tween a pair of antifermions (or a pair of each). The 
cr-exchange potential between quarks or quark- antiquark^ 



in the nonrelativistic problem in coordinate space can 
been written as 
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As for the 7r-exchange, since it has a different sign 
when taking non-relativistic limit from Dirac spinors to 
Pauli ones for the effective interactions u^uu^u and 
u^uv^v, so there is an additional minus sign from 
fcrmion-fcrmion into fermion-antifermion in taking the 
limit. 
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We note that the additional minus may be understood 
by the G-parity rule easily. Here ~<?i(j) is a Pauli matrix 
for spin and Tiu\ is a Pauli matrix for isospin. 

Since it is not allowed any quark-antiquark exchange 
between N and N, so the interaction from one-gluon- 
exchange for N and N is very different that for the NN 
system that one-gluon-exchange between N and N does 
not contribute to the effective potential at all. Here 
we only write down the potential due to one-gluon ex- 
change between two quarks (antiquark) in one baryon 
(antibaryon) [24] , 
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(ii). The effective potential derived from annihilation: 
Let us take it as an example, and for it the T-matrix 
of 7r-annihilation can be written as 

- "(pi, s' 1 )-i 5 v{p' 2 , s' 2 )-^gl h v(pi, s 2 )7 5 u(p,si) . (11) 

Using Fierz identities for Dirac matrices, a four- 
Fermion operator can be expressed definitely as a lin- 
ear superposition of others with a changed sequence of 
spinors as follows. 
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fc=l,16 



C ik {aO k d){cO k b) 



Since we are taking into account the contributions of 
the lowest order to the 'full' S'-wave at this step, i.e., we 
consider the 'full' S wave cases (all of the quarks and an- 
tiquarks are in S- wave), so under static approximation 
the contributions from cr-meson (J PC = ++ ) annihila- 
tion (P-wave annihilation) to the potential in the present 
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FIG. 3: The annihilation diagrams for quark and anti-quark. 



case (for S-wave states) are tiny. Thus we omit them 
here. The contributions from one-gluon annihilation to 
the potential between quark and anti-quark in momen- 
tum representation can be written as 

n = ^I( 16/9- l/3(Af • Af) 




1/2 + 1/2(7? -7?) - 3/2 + 1/2(5} • 5*) . (12) 



Here s = (pi +P2) 2 (in 4-dimension) in the propagator. 

In an annihilation process, s is time-like (s > 0), un- 
der static approximation so we have s ~ (mi +TO2) 2 , and 
when transferring the potential to space-time represen- 
tation, we have 
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4na' s 8(r) 1 



16/9 - l/3(Af • Xf) 



(rrii + rrij) 2 4 
1/2 + 1/2(7? ■ 7?) ) ( - 3/2 + 1/2(5$ • 5?) ) (.13) 



most parameters are fixed by fitting baryon spectrum 
of the model|2^|: The coupling constant g 2 h /^ is fixed 
b y 9nntt/ Ait , i.e., g 2 c Jkn is related to f% qir directly by 

^rZ^N &hl^ = fqqn> i m n,m N are the masses of pion 
and nucleon, respectively), the one-gluon-exchange cou- 
pling constant a s is determined by mass splitting of A 
and N, the confinement strength a c is fixed by the sta- 
bility conditions of TV, Vq is fixed by the masses of N. 
In summary, the parameters in the model are listed in 
Table I. The units for m q , m a , m n , Vq, b and a c are MeV, 
fm and MeV ■ fm~ 2 respectively. 

Tabic I. The parameters of the model. 



m q 


b 


a s a c 




m a 


m n 


A 


Vq 


313. 


0.60 


0.95 7.3 


0.59 


570. 


138. 


829. 


72.5 



To extend this model from NN systems to NN (nucleon- 
antinucleon) systems based on the effective Lagrangian 
Eq. © , now the Hamiltonians of the model may be writ- 
ten as 
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For 7r annihilation, we have 
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m 2 . - (rrii + 



1/2 - 1/2(5? • 0*) 



3/2 - 1/2(7? • 7?) 1/3 + 1/2$ • Af) (14.) 



where a c is the confinement strength. Tcm is the kinetic 
energy in the system of center mass. V qq elcch and V qq fxch 
are the effective potential from 7r and a exchanges be- 
tween the couple of quark and antiquark, respectively. 
V£ mn * and V* 
annihilation, respectively. 

Each of the wave functions for nucleon and antinucleon 
can be written as products of three parts respectively: 



and Vqq nn ' are those from the tt and one-gluon 



B. Hamiltonian and the wave functions of the 
model 

We follow the chiral quark model for multi-quark sys- 
tems [23], which essentially is an effective theory on ex- 
changes of Goldstone mesons, scalar meson a and gluons 
as well between quarks, but we extend it to antiquarks 
involved, so as to study the nucleon-antinucleon system 
with definite isospin (I) and spin (S). 

As the first step and to the conjecture of BES col- 
laboration here only S-wave states of the nucleon 
pair are considered, i.e., the total orbital angular mo- 
mentum L = 0, and we have J = S (the total an- 
gular momentum comes from quark spin only). The 
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For the orbital and color parts, the proton and antiproton 
have the same 'internal motion' wave function (here the 
motion for each center of mass has been removed): 



$0 = ( 2 N3/4( 2 \3/4 e -(Xi/(3b 2 )+pi/(ib 2 )) 
I 2 N3/4( 2 N3/4 -(X 2 2 /(3b 2 )+p 2 2 /(4b 2 )) 



$0 

P 



here A, p are the Jacobi coordinates of the components 
in each of the two clusters (nucleon, anti-nuclen), respec- 
tively. For the color factors: 

$p = -^=(ryb — rby + ybr — yrb + bry — byr), 

= "75 ^y^ ~ y^ ~ y ~^ y ~ ^ffi) ■ 
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For the flavor factors for a NN system, there are four 
possibilities with definite quantum numbers I and J. They 
are precisely (all symbols here have their usual mean- 
ings): 

\(jp T P I +P I P t -n T n i -n 1 n T) 

for I,J PC = l,l— ; 
|(P T P i ~P I P T -n T n i +n j n t), 

/or I,J PC * = 1,0-+; 
|(P T » 4 +P I P T +n t ^ I +n I n T), 

/or 7,J PC = 0,1— ; 

|(p T » 1 -p J. » T +n T ft I -n 4 ft T), 

/or /, J PC = 0, 0~+ . 

If we relate the states to the observation at BES with 
the decay J/tp — > pp+7[2^, we are sure that the C-parity 
of the pp pair must be positive i.e. C = +. Furthermore 
if we restrict ourselves to take S-wave into account only, 
then the states with minimal total angular momentum 
can be J PC = '"(Z = 0, S = 0) only. Whereas, here we 
also consider the states with J PC = 1~~(L = 0, S = 1) 
for comparison. 



C. A outline of the resonating group method 

For the the resonating group method [23|, first of all, 
to write down the two-cluster wave function with the 
conventional ansatz (to factorize out the relative motion 
of mass centers of the two 'clusters') as follows, 

|* rf > = [$p*p]M /5 x(2). (16) 

where [c] = [222] gives the total color symmetry. x(ii) is 
relative motion wave function of the two clusters. $ p and 
$p are the wave functions of the nucleon and antinucleon 
clusters in isospin and spin space only. 

To the specific problem, Gaussian functions with var- 
ious reference centers Si (i=l...n) are introduced, which 
(Si) play the 'generating coordinates' in the formalism, 

= (d^ )3/4e ^ { ~i (i?_5?i)2} 

and the relative motion wave function of the two clusters 
of the quarks and antiquarks is expanded into partial 
waves 

x(R) = "£x L (R)Y LM (k 

L 

N 

= ££ C fxf(^^ iM 0R), (17) 

L i=l 

with 

xf(R,s t ) = f dn s MR,&)Y LM &) 



where il is the L-th modified spherical Bessel function. 
For L=0, one has io(x) = sinh(x)/x (for a bound state). 

According to the ansatz of the RGM and having the 
center of mass motion 

<WJ? cm ) = (A)W(-^lp) 

7TO O l 

included, finally the wave function of six quarks within 
the two-clusters accordingly can be written as 

n „ 3 6 

k 1=1 J a=l/3=4 

here ip a (Si) and ippi—Si) are the single-particle orbital 
wave functions with different reference centers 

^(-^) = (^) 3 / 4 e-^+^/2) 2 /( 26 2 ). 

With the variational principle, one may obtain the RGM 
equation 

J H(R, R r )x{R')dR' = E J N (R R')x{R')dR' (19) 

via the variation with respect to the relative motion wave 
function x(i?). With a re- formulation, the RGM equation 
becomes an algebraic eigenvalue equation |27| 

j,k' j 

We should note here that for the nucleon-antinucleon sys- 
tem only the direct terms contribute. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

There are four possible states for an S'-wave nucleon- 
antinucleon system with different isospin / and total an- 
gular momentum J respectively. The spin, isospin and 
spin-isospin matrix elements of the interaction for the 
possible states are listed in Table II. 



Table II. The spin, spin-isospin and isospin coefficients 
for different states 



/ J PC 1 1— 


i o-+ 


l- 


o o-+ 


"1/9 


1/3 


-1/9 


1/3 


(n ■ fi) -1/9 


-1/9 


1/3 


1/3 


(<7i • ajTi ■ fj) 25/81 


-75/81 


-75/81 


25/9 
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For a qualitative analysis, firstly we neglect the 
nonlocal terms of the potential, so effective nucleon- 
antinucleon potential can been obtained by the Born- 
Oppenheimer approximation. The dependence of the 
'potential' (minus the 'rest masses' of the proton and 
the antiproton) on the separate distance Si is defined as 
the expectation value of the Hamiltonian at a separate 
distance Si as follows, 

(* PP -(SQ |ff 1*^(50) 

-{%\H\V P } - (y p \H\V p ) . (21) 

To see it precisely, we have calculated the nucleon- 
antinucleon effective potential with various gluon anni- 
hilation coupling constant. For simplifying, only the re- 
sults in the cases with I,J= 1,0 are drawn in FIG. 4. The 
contributions from 7r and cr-exchange, 7r-annihilation and 
one-gluon annihilation with different isospin and spin are 
shown in FIG. 5. 

In the calculations, all the model parameters are 
determined by fitting nucleon-nucleon interaction and 
deuteron properties (for the annihilation contribution, 
the parameters are related to those for exchanges through 
the Lagrangian Eq.lJBJ). From QCD, we know, one-gluon 
running coupling constant would become small with the 
increasing of momentum transfer. Here we let gluon an- 
nihilation coupling constant ot! a equal to one-gluon ex- 
change coupling constant as that a s in Table I. firstly, 
then reduce a' s appropriately to see the dependence on 
gluon annihilation coupling constants for the potential 
of the nucleon-antinucleon systems. For comparison, the 
case to ignore the contributions from annihilations are 
also computed. 

FIG. 4 contains 2 sub-figures. It shows the total 
nucleon-antinucleon effective potential with a' s — a s — 
0.95 and a' s = \a s = 0.316 for I, J = 1,0 state. FIG. 5 
contains 3 sub- figures. It shows the contributions from tt 
exchange, 7r annihilation, a exchange and one-gluon an- 
nihilation for the states I, J PC = 1,0"+; 0,1 — ; 0, 0"+ 
respectively. For the other states, we only present their 
minimum values of the effective 'potential' in Table III. 



Table III. List of the minimum of the effective potential 
for the possible states (unit for V m i n in MeV). 



IJ 


11 


10 


01 


00 


a£=0.947 


-20.1 


-18.1 


-14.1 


-14.9 


c/=0.4735 


-31.7 


-32.3 


-32.1 


-16.2 


c/=0.3156 


-37.4 


-39.7 


-43.6 


-16.7 


without annihilation 


-89.3 


-150.6 


-150.6 


-21.4 



If the attraction of the effective potential is deep 
enough for forming bound states, we further do dynam- 
ical calculations in the framework of resonating group 
method exactly, and finally obtain corresponding rela- 
tive motion wave functions and the mean square roots of 
the radius as well for the states. 
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FIG. 4: The effective potential for an S-wave nucleon- 
antinucleon system with different gluon annihilation coupling 
constant with state I,J = 1, 0. 

From FIG. 5 we may see that a exchange gives a uni- 
versally attractive interaction for a nucleon and antin- 
ucleon system as it does for a nucleon-nucleon system. 
As stated before, the contributions from a annihilation 
to the 5-wave nucleon and antinucleon system is small, 
so we omit it. ix annihilation gives repulsion interac- 
tion for all concerned cases with different strength, the 
contribution to I,J PC = 1,0 h is the largest, while 
the contribution to /, J PC = 0, h is the smallest. 
One-gluon annihilation also provides rather repulsion for 
these four considered states, especially, for the states 
/, J PC = 1, 0, 1 . (It seems that the contribution 
from the gluon annihilation should not be omitted when 
studying nucleon-antinucleon interaction, unless coupling 
constant of gluon annihilation is much smaller than that 
of gluon exchange.) tt exchange is the main factor to 
cause the differences for these four states. It provides 
repulsion for the state /, J PC = 0, h (and a minor 
attraction in the long range) but attraction for states 
/, J PC = 0, 1"; 0, 1 , that is opposite to the case for 
nucleon-nucleon systems. 

From Table III, one may see that the four possible 
states with different isospin and spin quantum numbers, 
all have intermediate attraction with different strength, 
and the attraction increases with decreasing the cou- 
pling constant a! s of the gluon annihilation. The ob- 
servation by BES collaboration presents a strong en- 
hancement in the decay J/tp — > jpp, that imply the 
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FIG. 5: The contributions to the effective potentiai from 7r 
and (T-exchange, ^-annihilation and one-gluon annihilation for 
the states I, J = 1, 0; 0, 1; 0, 0. Square points denote the con- 
tribution from gluon annihilation, circle ones for pion annihi- 
lation, up-triangle ones for pion exchange and down-triangle 
ones for a exchange. 



quantum numbers of the pp state, if it is in S-wave, 
are I (J PC ) = or 1 (0"+) but not I (J PC ) = or 
1 (1 ). When dropping the annihilation contributions, 
it is clear in Table III that the state with quantum num- 
bers I,J PC = 1,0 h has a larger attractive interaction 
than that with /, J PC = 1, 1"", that is consist with the 
BES collaboration observation. With the increasing of 
annihilation coupling constant, the difference of theses 
two cases will emerge. 



FIG. 6: Relative motion wave functions for the states 
I,J PC = 0,1""; 1,0"+ and 1,1"" with a' s = 0.3156 (the 
left figure) , and without annihilation contributions at all (the 
right figure). 



In summary, in this paper we extend to apply the 
chiral constituent quark model with quark-meson-gluon- 
degrees of the freedom to the nucleon-antinucleon sys- 
tems with quantum numbers (I,J PC ) — (1,1 ); 
(1,0"+); (0,1 — ); (0,0"+). Our results with the original 
model parameters as we can show the facts as below: 

(a) For the /, J PC = 0, h (S-wave) system, the re- 
pulsion from 7r-exchange is so strong that it can cancel 
the attraction from cr-exchange. No matter we adopt the 
annihilation coupling constant for the model in a reason- 
able range, there is no bound state at all. 

(b) The situation for the system with the quantum 
numbers I, J PC = 0,1"" and I, J PC = 1,0"+ for the 
,9- wave nucleon and antinucleon system. 

(c) By dynamical calculations we show that the bound 
states with /, J PC = 0, 1 ; 1, h ; 1, 1 for an S-wave 
nucleon and antinucleon system may be exist if the an- 
nihilation coupling constant a' s is suppressed for certain 
reason appropriately. For instance, if letting a' s — 0.3156 
for I,J PC = 1,0 ^, one may ob tain a 'tightly bound 
state' at E=1877.80 MeV and v r (r I ) = 2.06/m (see 
FIG. 6). If we had dropped the annihilation contribu- 
tions at all, except I,J PC = 0,0 ^ in the other three 
cases (I, J = 1,1 ;1,0 '"jO, 1 ) the nucleon and 
antinucleon system might be bound tightly within a size 
not greater than 1.1 fm. Moreover, the binding energy 
of the I, J PC = 0, 1"" and I, J PC = 1,0"+ states would 
be several tens of MeV greater than that of the state 
I, J PC = 1,1— and the mean squared root of the radius 
of the /, J PC = 1,1"" would be a little larger than that 
ofI,J PC = 1,0"+; 0,1"" etc. 

In the early 1990s, Dover et alfljjj constructed NN 
potential model from the NN effective potential by G- 
parity transformation accordingly, and predicted lower- 
lying isospin 1=0 natural-parity J PC = ++ , 1 , 2 ++ 
bound states and a few isospin 1=1 states, such as 
'"j 1 , also. Their predictions are consistent with 
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our results, if one compare those of an S-wave nucleon 
and antinucleon system only. It is interesting to note 
that there are so many similarities qualitatively, although 
Dover's model and ours are based on very different level, 
especially, to deal with the short-range behavior of NN 
in a very different way in the two approaches (in Dover 
approach, an arbitrary square-well cutoff is applied for 
the unknown short-range behavior of NN potential) . 

Finally, we should note here that this work is very 
preliminary in studying the nucleon-antinucleon interac- 
tion in the framework of constituent quark model. There 
are quite a lot of factors for the nucleon-antinucleon in- 
teraction which should be investigated carefully, espe- 
cially, as a very strong assumption, we ignore the couple 
channel (such as the multi-pion and other mesons be- 
ing involved etc) effects at all, although there are a lot 
of channels which should be considered |2^. Further- 
more, the adopted resonating group method also should 



be tested thoroughly. Even though, we still would like 
to emphasize that searching for the NN bound state 
with the quantum number /, J PC — 1, h through vari- 
ous multi-meson decay channels is crucial to confirm the 
multi-quark conjecture. Even though considering our re- 
sults and the BES observation, we would like to say that, 
if there is really an S-wave pp bound state, then its quan- 
tum number is likely to be /, J = 1,0 h . As a con- 
sequence, similar enhancements in the decays, such as 
B+ -> nnK+, B+ -> pnK s ; and B° -> nnD°, B° -> 
npD + etc in Belle and Babar at B-factories, and such 
as J/ip — > wry etc in BES at BEPC should be observed, 
although there are technical difficulties for the observa- 
tions. 
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